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SYNOPSIS 


Flow Induced Vibration is a major problem in pipelines, 
heat— exchangers arxi nuclear reactor cores , Reducing the level of th 
vibration has great technical importance in general, whereas in the 
nuclear reactor it gains its importance from the safety point of 
view also. To understand the basic mechanism behind this phenomenon 
considerable research work has been done in the past two decades. 
Still, some of the mechanisms are remaining fuzzy, especially 
that in the cross flow induced vibrations. 

Recently a mathematical model has been proposed for 
analysing cross flow induced vibrations in tube banks, using 
2-dimensional potential flow fluid forces. In the present work, 
this model is extended to the combined cross and axial internal 
flow. Using this model, the Finite Element solution for vibration 
of a tube bank in stationary outside fluid, as its special case, 
is obtained and it is comparable with the available earlier 
results. Results have also been obtained by FEM, for tube banks 
in (a) stationary fluid, (b) cross flow and (c) combined cross 
and axial internal flow, using the fluid forces obtained from 
potential flow theory and experiments, for a class of parameters 
such as outside fluid mass ratio, gap to radius ratio, number of 
tubes and axial flow velocity. In (b) and (c) cases, the potential 
flow theory predicts buckling mode of instability whereas the 
experimental results predicted dynamic instability. Use of 
experimentally measured coefficients also predicts the dynamic 
instability. So to establish close correlation with the experimenta 
results, the analytical formxalation for finding the fluid forces 
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CHAPTER -I 

INTRODUCTION 

1 .1 INTRODUCTION TO FLOW -INDUCED VIBRATIONS : 

The flow of fluid through and around the cylindrical 
structures, can cause destructive vibrations. These vibrations 
are induced by the fluid -structure interaction, i.e, the forces 
exerted by the fluid cause the structure to deform ; as the 
structure deforms, it changes its orientation to the flow and 
the fluid forces may change. This dynamic interaction may tend 
to increase in the vibration of the structure and can cause 
instability. 

Based on the configration of flow, the problems associate* 
with this phenomenon of vibration can be classified as having beei 
induced by 

(i) internal axial flow, 

(ii) external axial flow, 

(iii) cross flow and 

(iv) annular flow 

In case of axial internal flow, the fluid flowing throiigh 
the pipe induces the vibrations. Failures due to this vibration 
are reported in oil pipelines, he at -exchanger tubes and 
pump discharge lines [1,2 ], The mechanism of this flow induced 
vibration is well understood but the means of prediction of 
the excitation field are not yet adequate. 
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In case of external axial flow, the fluid, which induces 
the vibrations, is flowing outside the tubes and it is parallel 
to axis of the tubes. This type of flow occurs in heat -exchangers 
(parallel and counter flow) and nuclear reactor fuel bundles. 

The vibration characteristics of the external axial flow are in 
many ways similar to the axial Internal flow. Here also the 
mechanism of this flow-induced problem is well understood but 
the prediction techniques need further improvements , 

A large number of cross flow-induced vibrations are 
reported in transmission lines, tall chimney structures, closely 
packed tubes in heat -exchangers and nuclear reactor core [’1,2] , 
In these cases, the fluid flowing across the structures induces 
the vibration. Several excitation mechanism are possible in 
cross flow and they are 

1 . vortex shedding, 

2. fluid -elastic instability, 

3. turbulent buffeting and 

4. acoustic resonance. 

Except fluid -elastic instability, all other mechanism are 
possible in cross flow induced vibration of a single cylindrical 
structure. In this case, the mechanism of the vortex shedding 
with or without acoustic resonance is well understood but the 
turbulent buffeting mechanism needs further investigation. The 
available prediction methods are inadequate. 
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In case of cross flow-induced, vibra'tion in tube banks, 
several type of excitation mechanisms can occur simultaneously 
and it is difficult to identify which is dominating at \diat 
place. In fluid -elastic mechanism, the tubes whirl in an oval 
orbit at flow velocity above a critical velocity. The phenomenon 
behind this fluid-elastic instability is not fully understood. 
Even though some methods of prediction are available for this 
mechanism, they are inadequate, •^rom the failxxres reported 
under cross flow-induced vibration, [1 ], it is noted that most 
of the failures are caused by fluid -elastic instability. 

Some annular flov/ induced vibrations are reported in 
some major components of the nuclear reactor, e.g. the thermal 
shield and core barrel which lie concentrically within the 
pressure vessel, separated from one another by liquid annules [1 
A considerable amount of work has been done on the dynamics of 
cylinders containing or being immersed in annular flow. The 
fundamental mechanism for this vibration is fairly well urder- 
stood but the complexity of geome.-nry and flow fields make 
vibration prediction sometime difficult. 

In the present work, the effects of (a) cross flow 
(b) combined cross and axial internal flow induced vibration in 
tube banks have been studied, Sq, the previous research works 
which are related to this field are reviewed in the next 


section. 
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1 .2 REVIEW OF PREVIOUS WORKS : 

A comprehensive review about the research works done 
in axial flow-induced vibration are available in [2,3,4 ]. In 
early 1950, the vibration of a fluid conveying, pipe was studied 
in conduction with the axial flow-induced vibration of the 
trans-Arabian pipeline, '^he equation of motion of a fluid 
conveying pipe was derived by Housner [5 ]and he solved this 
equation for the trans -Arabian pipeline geometry . For the 
simply -supported pipe, he found that at low axial flow velocity, 
there was negligible effect upon the vibration of pipeline but 
at a certain high velocity, the pipe buckled like a column 
subjected to axial loading, Gregory and Paidussis [6 ] studied 
the axial flow-induced vibration in cantilever pipes and they 
showed that at sufficient high axial flov^^ velocity, the canti- 
lever pipes failed due to flutter type of instability. The 
effect of internal pressure of the fluid was studied theoretical 
and experimentally by Naguleswaran and Vvilliams [ 7 ]. In this 
case, they found that the effect was similar to the axial 
flow velocity, i.e, the pipe can become unstable at low -axial 
flow velocity by the action of interml pressure . Chen [ 8 ] 
studied the effect of displacement and rotational springs at the 
free end of the cantilever pipe. Using finite element method, 

Deb and Kumar [9f10 Janalysed this problem and studied the 
effect of different kinds of boundary conditions including 
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multiple supports on the stability of the pipes conveying 
fluid. 

The dynamics and stability of a pipe conveying fluid 
with pulsating components of axial flow velocity was studied by 
Paidoussis and Issid [3] , Mishra [11] studied the inpl^ine and 
outof plane vibration of cxirved and U - pipelines, using 
finite element method. Recently in 1986, and Kaladi [12] 
adopted the transfer matrix method for the study of vibration 
of pipelines , 

A brief review about cross flow -induced vibration is 
available in [1,13,1^] . In this, the vortex shedding mechanism 
has attracted the attention of investigators for many years , 

The vortex shedding behind a stationary cylinder can be 
described by the Strouhal number, S=2fg (^“ cylinder 

radius, U = cross flow velocity and f = vortex shedding 
frequency) . In case of a single cylinder, this number is equal 
to 0.2 and 0.3 for Re^ 2x10^ and Re > 3.5 x 10^ respectively [ 1 ]. 
(Re= Reynolds number) . In other ranges of Re , S is not sharply 
defined. An important phenomenon associated with vortex shedding 
excitation is lock-in. i.e., the vortex shedding frequency 
may concide with the natural frequency of the cylinder. In this 
lock-in region, the vortex shedding mechanism causes the 
instability . 

Earlier researchers had doubts about the existence of 
vortex shedding in closely packed tube banks, because of 
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insufficient space for the shear layers to roll up in classical 
manner. But later investigations proved that the vortex 
shedding was possible in atleast first few rows [15] and 
Strouhal numbers associated with flow across different layout of 
tube banks are available [ 2] . Pettigrew and Sylvestre [15] 
studied this vortex shedding ejfcitation in liquid cross flow 
tube arrays , This phenomenon appears to be absent in two phase 
flow and it remains fuzzy . There is a possibility of the 
concidence of vortex shedding frequency with acoustic of the 
flow channel and the natural frequency of the tube banks . Such 
a case may lead to the damage of the he at -exchangers [ 1], 

Connors [16 ] was the first investigator to recognize 

the fluid -elastic mechanism of vibrations caused in the tube 

banks subjected to cross flow. Based on the experimental mode 

shape, he measured the quasi-steady fluid -elastic forces and 

developed a simple stability criteria for a tube-row subjected 

U 

to cross flow in terms of reduced cross flow' velocity, 11^,= ■ > 

and mass -damping parameter, d = , (f== natural frequency 

of the tube, t = damping ratio, m= mass per unit length of the 
tube and p = density of the cross flow fluid ) . 

U * k d (1.1) 

r m 

where k is a dimensionless instability constant and is a 
function of fluid forces which in turn depend on array geometry. 
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He obtained the quasi -steady fluid forces by assuming that 
the dynamic fluid forces acting on a tube at any instant were 
given by steady fluid forces acting on that particular 
instantteous position and its neighboiiring tubes. Connors [16] 
found that k « 9.9 for his tube-row and Eqn. (1 ,1) has been 
extensively used as a design guideline for the prevention of 
fluid -elastic instability of heat -exchanger tube array, 

Belvins [ 17 ] extended the model to include the variable 
tube-to-tube stiffness and damping. He also attempted to find 
an expression for the stability constant of a full tube array, 
using theoretically determined fluid forces [18 ] , In his book 
[2 } ,he suggested that the additional flow velocity -dependent 
damping terms should be included in the equation of motion of 
the tubes . 

In the above models, the stability curve was found by 
assumed mode shape of vibration of the tubes. Price and 
Paidoussis [19] have relaxed the assumption of mode shape, 
instead the measured fluid force coefficients are inserted into 
the equations of motion of two-row array and stability is 
examined using an eigen value analysis. In this way, the least 
stable mode of vibration is predicted, it being the eigen vector 
corresponding to eigen values of the lowest critical cross flow 

velocity . 
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The mathematical approach iised by [16,18,19 ] Includes 
multiple tubes with fluid -elastic and fluid -damping couplings 
but the fluid -inertia coupling, which is important for closely 
packed tube banks vibrating in dense fluid, is not accounted. 

Chen [ 20-22] studied the fluid -inertia coupling effects 
on vibration of group of circular cylinders in [a) unboionded 
and (b) bounded ideal still fluid. Using 2-dimenslonal 
potential flow theory, he found the fluid -inertia coupling 
between the tubes. He extended the potential flow theory to the 
cross flow problem, [23], for finding out the fluid -inertia, 
fluid -damping and fluid -elastic coefficients . These coefficients 
are dependent upon the array geometry and are independent of 
cross flow velocity. 

Tanaka and Takahara [ 2A] were able to experimentally 
determine these fluid force coefficients for water flow across 
a row of tubes and a square array and found them to be dependent 
on the cross flow velocity as well as the motion of the tube 
in question and its nearest neighbours. They were able to 
perform a nximerical stability analysis to determine the least 
stable mode of vibration and generated a stability curve covering 
the region of their water test. Chen [25] used the Tanaka -Takahara 
force coefficients and extended the stability curve into the 
air -flow region. Some interesting jumps or discontinuities in 
the stability curves were observed by both group of researchers. 
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Number of conferences have been held in the flow induced 
vibrations area and the papers presented in some of the 
conferences are available in [28-31 ], 

1 .3 OBJECTIVES AND SCOPE OF PRESENT WORK : 

The present work describes the finite element analysis 
of tube banks vibrating in (a) stationary outside fluid, 

(b) cross flow and (c) combined cross and axial internal flow. 

The finite element formulation of the flow induced vibration 
uses the standard Galerkin finite element method. 

The finite element method is tested with the vibration 
of simply -supported tube-row in a stationary fluid, using 
2-dimensional potential flow fluid forces. The finite element 
solution is compared with that of the closed form solution of 
Chen [20], Also the effect of parameters (a) outside fluid mass 
ratio, V, (b) gap to radius ratio, G/R, and (c) the number of 
tubes, K on the natural frequencies of the tube hanks are studied. 
In addition to this, different type of single span and multispan 
tube -rows are analysed. 

In the case of cross flow, the critical cross flow 
velocity at which the stability occurs first is studied for 
different layout of tube banks , The fluid forces acting on the 
tubes due to cross flow are obtained from 2-dimensional potential 
flow theory. Experimentally measured fluid forces are also used for 
analysing tube-row. The stability of the system is determined 
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by the eigen values of the system. Lastly the combined cross 
and axial internal flow is studied for a row of 5 tubes. The 
potential flow and experimentally measured fluid forces are 
used for tlie stability analysis. 

In Chapter II, the equations of motion for the combined 
cross and axial internal flow are derived and its finite 
element formulation is given. Formulation for cases cited above, 
i.e., (a) stationary outside fluid and (b) cross flow are 
special cases of this formulation. 

Chapter-Ill presents results for all the three cases 
mentioned above. For the stationary outside fluid, parameters 
considered are outside fluid mass ratio, gap to radius ratio 
and the number of tubes in a row. For the cross flow^ e feet of 
mass -damping parameter on the critical cross flow velocity 
has been studied. Whereas for the combined cross and axial 
internal flow, the critical cross flow velocity has been found 
out for different mass -damping parameter and axial internal 
flow velocity. Based on the results and discussion, conclusions 
are derived and are stated in Chapter IV . 
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CHAPTER II 

EQUATION OF MOTION AND FEM FORMULATION 


In this chapter, the equations of motion for a 
bundle of circular tubes subjected to the combined cross flow 
and axial internal flow are derived. To obtain the solution 
of these simultaneous partial differential equations for 
various boundary conditions and multiple supported tubes, 
Galerkin finite element method is used to convert these into 
simultaneous oixLinary differential equations. 


2,1 EQUATIONS OF MOTION : 

A bundle of K tubes subjected to cross flow is shown 
in Fig,1. The axes of the tubes are parallel to Z- axis and the 
cross flow is parallel to the ^-axis with velocity U^ . The 
subscript i is used to denote variables associated with tube 
i. Thus E. I^ and are the flexural rigidity and mass 
per unit length of the tube i. 

The equation of motion for tube i in the X direction is 


k 

d u 


^i^i” 


. ®^^i 

+ <=1 w 


i » 1,2,3, ...,K (2.1) 


where u^, and f^ are displacement, structural damping 
coefficient and fluid forces per unit length in x direction 
Similarly, the equation of motion in Y direction is 
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Fig.l Bundle of Tubes in Cross Flow. 
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av. 3^. 

TT" ”Si » i-1,2,3, ...,K (2.2) 

where v^, and are displacement, structural damping 
coefficient and fluid forces per \init length in Y direction. 
With the following assumptions ; 

1. The fluid flowing across the tube bank is inviscous, 
incompressible and unbounded. 

2. The cross flow is irrotational one, 

3. The amplitude of tube motion below the stability 
threshold is vei^ small. 

4. There is no separation of flow. 


the fluid forces f. and g. are calculated for tube 

i i 

banks arranged in arbitrary pattern, using two dimensional 
potential flow theory, Chen [23, 25 ] . In this approach the flui< 

velocity potential is expressed as a truncated series of 
unknown coefficients. The boundary conditions are imposed by 
transforming, the velocity potential to the tubes natxiral 
coordinates. The unknown coefficients are computed by solving 
set of equations which are obtained by imposing boundary 
conditions. The fluid forces are calculated from fluid pressure 
acting on the tubes and are. 

K R.+R. 2 3^u. 3^. K ^ 3u. 3v. 

V (“Id 
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+ pU 


o 


K 
2 








+ PR, U‘ P, 


(2.3) 


>#here 

“i;)* 



ij* io* i;3 


added mass coefficients, 
fluid damping coefficients, 

fluid -elastic stiffness coefficients, 


- X component of steady state fluid foroe coefficients, 

- Y component of steady state fluid force coefficients , 

p - density of cross flow fluid, 

U - cross flow fluid velocity, 

o 

R - outside radius of tute. 

The equations of motion for flexural vihration of a uniform 
tube i subjected to axial internal flow are, [ 3 ], 




3^u> ^^^i 2 . ^^^i 




aV 2 ®^i 


0 


0 


(2.4) 


where 


“ai 


mass 1 ) 6 r unit length of internal fluid. 
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From Ekins.(2,1) to (2.4), one can get the equations of motion 
for a bundle of tubes sub;)ected to the combined cross flow u , 
and axial internal flow V . , as follov/s 




Zl 


8Z 

+2m„ .V . 


Tf 


+m„ . V 


2 


K 


ai\i-T:3 -PK S (a- . u,-.a'' v,) + c, 


3u. 


ai ai ^TTE j ^ 3^) + (V“ai^ ^ 


32- - j=1 5 'O' 


i Ft 


, 2 , , 

c u. 


K • R,+R^ 2 


ix *v * -rjit . Cm O U» 

tjl|3 S ( 5J— U) (q; ^ + 


0=1 


id "■'^id 7 ^ 


i^)-pR.U" = 0 


j2 

'■±'^0 "^i 


at' 

a^^ K 3u . av . a^, 


(2.5a) 


3 . - 3 ^ 

1 . _ i ^Tr2L 


9 V. K au, 3 v. a'Sr. 

+PRi°o 7^ *^id 

K R.+R. 2 a\i. a^. 

+KP E ( — 5~*i) ^ _.^)-pR^u2 ^ Q (2.5b) 


d=i 


id 3^2: '“^id 


i = 1,2,3 ,K 


In the present work, all tubes have been taken as 
uniform and identical. Thus length, radius, mass, flexural 
rigidity and structural damping coefficients of all tubes are 
equal . Also fluid flowing through and its velocity are same 
for all tubes . So 

\ ' “i = “ * “al " “a 




a’ 


EiIj^=EI; 




(2.6) 
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2.2 FEM formulation : 

All tubes are divided into (q-1 ) elements with q nodes, 
see Fig. 2. Typical finite element is also shown in the same 
figure. The displacements u^*^and over the typical element 

of a representative tube i are taken as 

/ \ a ( ) 

u}®^(Z,t) = E N^(Z) u^(t) = LNjCUj_(t)l 

r«1 (2.7) 

v^®^(Z,t) = E^Mj,(Z) v^(t) =[MjW.(t)3 


where 

a - number of degrees of freedom associated with u displacemen 

b — number of degree^ of freedom associated with v displacemen 

[N J - shape ftinction for u displacement 

[M J - shape function for v displacement 

Substituting Eqns. (2.7) in Eqns . (2.5) residues Ru^ 
and Rv^ are obtained. These residues are minimized by the 
Galerkin method, thus 


h 

/ Ru^ dZ 


h 

/ Mp Rv^ dZ 


0 


0 


i “ 1>2, 


( 2 . 8 ) 


Eqns, (2.8) are integrated by parts as much as possible. 
It is seen that the highest order of derivatives of u and v 
are two in the integrals and three in the whole equations . 
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Nodes and Finite Elements of Tube i 



Typical Finite Element 
Fig. 2 Finite Element Representation- 
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-To ensure monotonio convergence, the element should 
be used for both u and v . The shape function for this 
elements are given as , 


[NJ= LMJ= L J(2-3^ ^(2+3^ 

|( j 


where ^ 


(Z^- Z )/2 


h = z^-z^ 


( 2 . 9 ) 


Using these shape functions, Eqns, ( 2 . 8 ) which have 
been integrated by parts can be written in matrix form as 

h \ o ^ 

El / ? N*J [ n'j dZ C u 3 ^ ^ -®a^a / ^ N jdZ C 

o ^ 

-PUq ^ (c^ij / {N}[NJ dZ {u^ / ?N}[Nj dZ ^ ^ ^ 

0=1 o ° 

+c ; C N) [N J dZ { ^ + 2 mg^\ { N} [N J dZ { u^} 


K h (ne) ^ , f * a V 

+ pRU^ r (a^^ / {N}LNjdZ{u^ +o'^. / {N}lNJ dZ I v J ) 

h (ne) 

.(m.m,); lN)lNJdZCU,)^-) nUN J dz Tu .i 

h (ii® ) 9 ^ , 

+ 0 . . / US lN J dz { v^} )- pRU^ / C N ) dZ 

JL^ 


+ {N} El 


3 „(e) 


a"u 


SZ' 


- i 


3 ^j 


El — K 

az^ 


(e) 

+{nJ “a|- 


au 


J oi 


( 2 . 10 a) 
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^ » h 

El J Cn) N J dZ{v.i [N | dZ{ v, } 

« 1 a a ^ ^ i 


K 


-■ PUq ^ (•^ij I ^N}iNJdZ{u^}'^^®^ + ^"^ /{ H|NjdZ {v^}'^^®-') 

+c ;{ nHnJ dZ{v^l^^®^ + /{ NJ ',N J dZ 


(ne)_ 


,(ne) 


K 


+pRU^ I ; {N} [NJ dzcu^j^’^^^ + p^^ /{ N}LNJdZ ) 

o o 

+ (m+m ) ; { N} [Nj dZ{v. 3^^®^+npR^ 1(t .. ; !N} [N J dZ { u . } 

o ^ ;3=1 o ^ 

li h 

+ ^. . ji NJ [Nj dZ! -pRU^ P. ;{ N} dZ 

J o ^ ^ ^ o 


aV®> 

+ {NJ ei ^ 

8 Z^ 


h ^ gS^Ce) 

-IN }EI i- 

o az^ 


av^®) 

2 


+ jN}HiaVa 32 — 


^ 0 ^ 


(2.10b) 


Eqn, (2,10) can be written as 

+®id["cR] )- [Fgjj) 4 !B^ 


(2.11a) 
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K 

Z 


Je) 


[KcRj'" [KaI 


i(®>- 


(e) 


(e) , f • 1 (®) .. 

■"^idE'^Al -"PidLCcRl ) +’'ld[«CR] tUjl 


, (e) (e) (e) (ne; 

[^ 1 •*•^15 [«a] "^i;3^CRl 


(e) (e) 


(2,11b) 


where 6 = Kxonecker’s delta 

(e) h . „ 

[ICj] = El j C N } [N JdZ, stiffness matrix due to tube bending ; 
o 


(e ) P h , , 

[K. ] = / CNiLNjdZ, stiffness matrix due to axial 

o internal flow ; 


(e) 


[K^^r] “'PUq / CN}[NJ dZ, stiffness matrix due to cross flow ; 


(e) h 

[ ] = c / { N } i^N j dZ , damping matrix of tube ; 

o 


r ^ r - ^ 

[C, J “2m V / tls&lNJdZ , damping matrix due to axial 
^ ^ ® o internal flow ; 


(e) h 

[^CrI i damping matrix due to cross flow ; 


o 


(e) h 

j^Hj] = a / {Ni[Nj dZ, mass matrix of the tube ; 


(e) h 

[M^] " ®a / {NiLNj dZ, mass matrix due to axial internal flow ; 


(e) « b 


[^r] / C N} lNj dZ, mass matrix due to cross flow ; 

o 
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(e) 


h 


P j 

Q^N}dZ, force matrix due to cross flow ; 

^ 3 ^ 

a^u. . 

El — 5^ e. 

a a 1 


'®u> 


(a) 




-EI- 




37' 
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3 u. 5 

-^^—3 -“a^a 
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3^1. 
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^ 3Z- 

^ aV 


,X-direction forces at the nodes. 
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^ Va «'i 








-EI- 


3Z' 

3 

8-^v 








-EI— ^ -m i. 

8Z ® ^ 

3^. 

El ' X- 

3Z^ 


,Y-direction forces at the nodes. 

( 2 . 12 ) 


On assembling, Eqns . (2,11), equations of motion for 
the i-th tube become 

I ( [0U1] [0V1] ^jCv.! 


+ [MU1 ^ ‘^3 ^ lij ^ ^ =^ FU } ^ + { BU } 
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[KU2]ij tUj! ^ [KV2] . [CV2] iv^! 

+ [m2-]^,{\l.) +[ HV2] CVjl ) = {FV} (2.13) 


where 
[KLJ1 ] 


t- -> i;3 “ ^ ^ ' 

(e) (e) „ (e) 

[KV2],^ = E(6.. [K^] -Pi^[Kcf,] ) 

, (e) (e: 

rmri 1 „ v-n* rv i . ^■[/^To i _ v_^' i 


, (e) (e) 

^^"^31^ “ (^^r 3 » 3 i j= [^r3 

(e) (e) (e) 

[CUI 3 i [C^] + 6 ^^ [C^ ] + 0 ^^^ 

'-) (e) . (e) 

■^*^i;3{^^A3 L^cr3 3: 


[CV2]^^= E(6.^[C^/ "'=±0 L^CR 

[CV1] , .=20^ . [C,^ l"\ [CU2] , .=2 T [C^^ i"\ 

(e) (e) (e) 

[MU1 = 2(6^^ ] +6^^ [M^ ] ] ) ; 

(e) (e) (e) 

[MV2],. =rC6i^[M^] [MchI ) # 


rmri]i^=2 i[MU2]i^ = 

(e) 




(e) 


FU} ^ = E 

(e) 

C{BU ) 

i 


BUK = I 


; CFV}^= EB. IF^J,} 

(e) 


(e) 


(2.14) 


; { BV } = E{ BY } 

i i 

As said above the Eqns . (2.13) are for a typical tube i 
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where 
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( 2 . 16 ) 


(2.17) 


Eqn. (2,16) consists of simultaneous, second order differential 
equation which can be solved for various type of boundary 
conditions . 
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CHAPTER -III 

RESULTS AND DISCUSSIONS 

In this chapter, the method of solution for the set of 
differential equations, which are obtained from equations of 
motion by FEM, is discussed. Then the results for the following 
three cases are presented, 

(i) Vibration of tube banks in stationary fluid 

(ii) Vibration of tube banks in cross flow 

(iii) Vibration of tube banks subjected to the combined 
cross and axial internal flow. 

The effect of different parameters on these results are 
also discussed, 

3.1 METHOD OF SOLUTION : 

The solution for the simultaneous second order ordinary 
differential equations, Eqn, (2,15), is obtained as explained 
in Meirovitch [ 32 ] and Duncan and Collar [ 33 ] , Using generalised 
velocity of and {v J as auxiliary variable, n second order 

differential equations, Eqn. (2.15), are converted to a set of 
2n first order differential equations . 

+[kJwJ = Cf > (3.1) 

where 







Cu} 

tvl 
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For damped free vibration, eigen values are obtained by 
assuming the following solution 


C W i = e^^ I W J 


(J.2) 


where K is an eigen value and is the corresponding eigen 

vector. Substituting Eqn,(3,2) in the homogenous form of 
Eqn, (3,1) leads to 

X riil^ wi + [K]^w5 = {03 (3.3) 

or 

[D]f v7l = 1 [I]SW3 


Where 


[ D> -[K] 


[0] [I] 


The eigen values of the dynamic matrix [D] , will in general be 
a complex number. The real part of the eigen value determines 
stability of the system and imaginary part gives the natural 
frequency of vibration. If X is complex, having negative real 
part, the system is stable and performs damped oscillations. 

When X is complex, having positive real part » the system becomes 
xinstable dynamically (flutter). When a is real, then the system 
loses stability by divergence (buckling). 



26 


3,2 DIMENSIONLESS PARAMETERS : 


The results are presented in terms of dimensionless 
parameters used in literature and they are 


Inside fluid mass ratio 
Outside fluid mass ratio 
Mass -damping parameter 

Dimensionless Axial internal flow velocity 
Reduced cross flow velocity 


P = 

V = 

6 = 
m 


m+m„ 

m+m„ 

tit ) 


V = 


U_ 


V^L 

a 


“o 

zTr 


Real part of dimensionless frequency 


= Re(D) = 


^El/((m+in^)L^) 


Imaginary part of dimensionless frequency 


=Im(D) 


-Sn(\) 

jEl/((m+m^)L^) 


where 


L - length of tube, 

f - fundamental natural frequency of the tube in 
vacuum (Hz), 

T - damping ratio (t = 4^ (m+mj ' f > 


3.3 NUMBER OF FINITE ELEMENTS USED IN EACH TUBE ; 


In finite element method, the accxiracy of the natural 
frequencies depend upon the number of finite elements used in 
each tube. To decide the minimum number of elements required for 
each tube, natural frequencies of a row of 5 simply-supported tubes 



27 


are obtained for 1 to 6 number of elements in each tube. 

The first and. second natural frequencies are given in Table 1 , 

It was found that after 4 elements, there was hait’ any change 
in natural frequencies corresponding to first mode of single 
tube in vacuum, S© all results are obtained by taking 4 elements 
in each tube, 

3.4 VIBRATION OF TUBE BANKS IN STATIONARY FLUID : 

This case is studied without axial internal flow and 
damping. Equations of motion are obtained from Eqn, (2.15) by 
substituting = m^^ =» = c = 0. It is noted that there is a 

coupling between the tubes in the mass matrix, [m], and it is 
due to the added mass coefficients . These coefficients are 
obtainable from the Z-dimensional potential flow theory. 

First, the natural frequencies for a number of simply 
supported tubes (1 to 10) in a row are obtained, gap to radius 
ratio, G/R, and outside fluid mass ratio, v , are taken as 0,2 and 
0,54267 respectively. These results are shown in Fig, 3 as a 
function of number of tubes in a row, K, It is known that 
corresponding to each mode of a single tube in vacuum, there 
are 2K coupled frequencies for a row of K tubes, [ 21 ] . Results 

in Fig. 3 are first, second, 2K-1^ and 2K^ frequencies corres- 
ponding to the first mode of single tube in vacuum. The mode 
shapes for these frequencies are also shown in the same figure , 
With increase in number of tubes, the first and second frequencies 
decrease and 2K*-1^^ and 2K^^ frequencies increase. All other 
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Table 1 : NATURAL FREQUENCIES OF A ROW OF 5 SIMPLY SUPPORTED 

TUBES IN STATIONARY FLUID FOR DIFFERENT NUMBER OF FINITE 
ELEMENTS IN EACH TUBE (G/R = 0.2, a-^0.5, 0,0) 


Number of 
elements in 
each tube 


1 

2 


3 

4 

5 

6 


First "5'requencv 

■ReCn) -T---^Frj27 


Second ^'re 

TTiPi r 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


7.0419 

6.3696 

6.3497 

6.3462 

6.3452 

6.3449 


0.0 


0.0 

0.0 

0.0 

0.0 

0.0 


7.7067 

6.9708 

6.9491 

6.9452 
6 .9442 
6.9^38 






Fig. 3 Effect of Number of Tubes on Natural Frequ^cy of a 
Row of Tubes in Stationary Fluid . 






5o 


frequencies corresponding to first mode of a single tube lie 
between the first and 2 k"^^ frequencies. Results for the first 
and second frequencies are given by Chen [ 20 ] and the present 
results match well with those results. 

Next, the effect of fluid mass ratio, and boundary 
conditions are studied and the result’s for the first natural 
frequency are given in Fig, 4, For mass ratio study, all tubes 
are taken as simply supported. With increase in mass ratio, 
all frequencies decrease, as expected. For boundary condition 
effects, all tubes are taken as fixed -fixed. All frequencies 
increase considerably as expected. Fig, 4 shows, these effects 
on the first natural frequency. 

Natural frequencies also depend upon the gap to radius 
ratio, G/R. This effect is studied for a triangular array of 
3 simply supported tubes . The arrangements of the tube array 
and results for the first six frequencies corresponding to first 
mode of single tube in vacuum are shown in Fig, 5. The gap to 
radius ratio, G/R» is varied from 0.01 to 1.0. As the G/R ratio 
increases the first three frequencies increase, the effect being 
very much pronounced in the first frequency. For the other three 
frequencies, the trend is opposite. As the spacing increases all 
the frequencies approach the first mode frequency of a single 
tube vibrating in infinite fluid, which is shown in dotted line 
in the Fig, 5, These results also match well with those of 
Chen [ 21 1 .The same trend is also observed in the set of 
six frequencies corresponding to higher modes of single tube in 
vacuum, though they are not presented here. 
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Lastly the effect of intermediate support on the 
natural frequencies is studied. The first two natural frequencies 
of hinged -hinged, fixed -fixed and fixed^inged tubes, 5 in a 
row with and without intermediate support are given in Table 2. 
Introduction of intermediate support in a tube row increases 
the frequencies, as expected. For these three boundary conditions, 
the effect is more pronounced in the hinged-hinged tube row, 

™MTION OF TUBE BANKS IN CROSS FLOW : 

The equations of motion for tube banks vibrating in 
cross flow can be obtained from Eqn. (2.15) by substituting m 

^a zero. There is a coupling between the tubes in [ M 1 ,[C ] 
and [It] matrices due to added-mass, fluid -damping and fluid- 
elastic coefficients , The stability of the tube banks is determined 
by the eigen values of the system. The instability can be static 
(buckling) or dynamic (flutter). The cross flow velocity, at 
which the first instability occurs, is noted as critical cross flow 
velocity , 

In this stability analysis, the damping ratio, t, and 
outside fluid mass ratio, Vt can be taken as two independent 
parameters , But Chen [ 25 ] and Blevins [ 2 ] found that both 
parameters appears in the stability criteria as a single parameter, 
mass -damping parameter, ^hey presented the results in terms 
of this parameter. Same has been followed in this work. 





TABLE 2.: NATURAL FREQUENCY OF A ROW OF 5 TUBES IN 
STATIONARY FLUID FOR DIFFERI'.N’T BOUNDARY 
CONDITIONS, j>= 0.5, G/R » 0.2. 
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First, stability of a row of five uniformly spaced, 
simply- supported tubes is studied. The gap to radius ratio, G/R, 
and the damping ratio, t , are taken as 0.66 and 0,02 respectively. 
The arrangements of the tubes and flow direction are shown in Fig ,6. 
By 2-d.imensional potential flow theory, the fluid forces acting 
on the tubes are calculated and are given in Table 3 to 5. 

Using these fluid forces, the critical reduced cross flow 
velocity, , at different mass -damping parameter, , are 
computed and plotted in ^ig, 6. It is noted that the tubes 
become unstable by buckling and the critical reduced cross flow 
velocity, , increases with 6^. 

Stability of a square array of 4 tubes and a triangular 
array of 3 tubes, simply supported, are also studied. The results 
for these cases are presented in ^ig. 7 and 8. 

Tanaka and Takahara [24 ] experimentally measured the 
fluid forces coefficients for a row of tubes subjected to cross 
flow of water. These measured coefficients are not matching with 
the potential flow solution and also noted that these coefficients 
are dependent of cross flow velocity whereas those 
obtained from the potential flow theory are independent of 

cross flow velocity . 

Chen [ aS] studied the published experimental data of 
critical cross flow velocity for various array of tubes . He 
observed that tubes were becoming unstable by flutter and not 
by buokllng. Based on this oolleoted Informations. he gave the 
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Table 3. 

ArifipH -mass Coefficients for a row 

of 5 tubes 


based 

on potential 

flow solution 

• 



1 .113 

-0.282 

-0.282 

-0.04 

—0 .04 


-0.282 

1 ,112 

-0.037 

-0.287 

-0.017 

[ “ia b 

-0.282 

—0 ,037 

1 ,112 

-0.017 

-0.287 

—0 .0^ 

-0.287 

-0.017 

1 .058 

-0.011 


-0,04 

-0.017 

-0.287 

-0.011 

1 .058 



■= [0 1 





1.135 

0.34 

0.34 

0,121 

0.121 


0.34 

1.13 

0.13 

0.319 

0,064 


0.34 

0.13 

1 .13 

0.064 

0.319 

0.121 

0.319 

0.064 

1 .065 

0.035 


0.121 

0.064 

0.319 

0.035 

1 .065 
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Table 4, : Fluid. -damping coefficients for a row of 5 tubes 
based on potential flow solution. 




[ ] 


0.0 
-1 .266 
1 .266 

-0.344 
0 .344 


1 .221 

-1 .221 

0.244 

-0.244 

-0.174 

—0 . 28 

1 ,078 

-0.107 

0.28 

0.174 

0.107 

-1 .078 

-1 .354 

-0.155 

-0.909 

-0.078 

0.155 

1.354 

0.078 

0.909 


r 1 


0.0 
—0 .93 
0.93 
-0.003 
0.003 


0.974 

-0.974 

0.04 

-0.134 

-0.032 

0.98 

0.032 

0.134 

0.008 

-0.741 

0.006 

—0 . 426 

-0.006 

0.741 

-0.005 


-0.04 

-0.008 

- 0.98 

0.005 

0.426 
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Table 5.: Fluid -slastic coefficients for a row of 
5 tubes based on potential flow solution, 


' ' -1 


5.752 

-2.6 

—2 .59 

-0.281 

-0, 

.281 

-2.591 

5.356 

—0 .265 

-2.4 

-0, 

.098 

-2.6 

— 0 ,266 

5.354 

-0.098 

-2 

.391 

-0.28 

- 2.393 

-0.098 

-2.84 

-0 

.06 

-0.281 

-0.098 

-2.401 

••0 #06 

2 

.831 


r 0 ". 1 = r x!'. ] = [0 1 
L 10 J '13 ^ - 


f bp 


•3.03 

1 .507 

1 .507 

0.007 

0.007 

1 .508 

-2.802 

0.014 

1 .278 

0.0004 

1 .507 

0.014 

-2.801 

0 .0004 

00 

CM 

0.007 

1 .286 

0,0004 

-1 .296 

0.001 

0.007 

0 .0004 

1 .286 

0.001 

- 1 .295 



Critical Reduced Cross Flow Velocity, Ur. 
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5 6 8 10 20 30 4050 t0( 


Mass - Damping Parameter, 6m 

Rg. 6 Critical Cross Flow Velocity ut, for a Row of 5 Tubes 
in Cross Flow. 


Critical Reduced Cross Flow Velocity, Ur- 
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Mass -Damping Parameter, 6m 

Fig. 7 Critical Cross Flow Velocity uT.for a Square Array 
in Cross Flow . 


Critrcal Reduced Cross Row Velocity, 
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Mass-Damping Parameter, 6m 

Critical Cross Flow Velocity uT, for a Triangular Array 
in Cross Flow. 
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empirical correlations for the lower bound of critical reduced 
cross flow velocity, [ 26 ] . Using these correlations, the lower 
bound of critical reduced cross flow velocity for the above 
cases are plotted in Fig. 6 to 8. 

Using the measured fluid force coefficients of Tanaka, 

Chen " 25 ] studied the stability of a tube row. It was 
observed that there was dynamic instability of tubes. This can 
arise due to fluid damping or fluid-elastic forces. In case of 
fluid dampin£- forces, energy transfer from cross flow to tubes is 
by tube velocities and the instability so caused is attributed 
to velocity mechanism and fo.r fluid— elastic forces, energy 
transfer is by tube displacements and the instability so caused 
is attributed to displacement mechanism. 

Vdien cross flow velocity is small, the velocity of the 
tubes will be of the same order of the magnitude of the cross 
flow velocity, velocity mechanism causes the instability .Chen r25j 
found that this happens for U^< 15 for a row of 5 tubes and 
experimental results showed that the tube oscillation are 
predominant in lift direction. 

When cross flow velocity is large, fluid-elastic forces 
are dominant and displacement mechanism causes the instability. 
Chen [25] found that this instability can occur for U^, > 15 

. c For the row of 5 tubes, his experimental 

for a row of 5 tubes, -"or -cne low u 
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results showed tliat the middle tubes oscillated predominately 
in tlio drai’; direction while the neighbouring tubes oscillated 
predoininantely in the lift direction. 


The critical reduced cross flow velocities, U^, » fo^ ® 
row of five tubes are calculated by FEM, using the measured 
coefficients of Tanaka [ z4] and model of Chen [25]. ^or 
tubes were taken as vibrating in lift direction, and for 
middle tube was taken as vibration in drag direction, next two 
in lift direction and last two as rigid, Chen [ 25 ] .He also 
neglected the fluid -damping force, fluid interia force and 
diagonal el'emunts of the fluid -elastic force for > 15. These 
results are shown in Fig. 6 and match well with those given ty 
%cn r 1 , It may be noted that if flxiid -damping, fluid inertia 

and diagonal elements of the fluid-elastic forces were n 
neglected, the tube-row was found to be stable. For this 


5 tubes, it is seen that dynamic instability is caused by 
velocity meohunism for 6^ < A and ty displacement mechanism for 

The critical reduced cross flow velocities, , obtained 
by potential flow fluid forces and measured fluid forces di^f ^ 
widely. The mode of failure is buckling based on potenti 
theory and dynamic instability based on measured fluid force . 

It can be explained as follows . 


44 


Imtubillty can ariae due to lluld-alastlc or fluid 
da,„pln,. forcaa. Table 5 gives the potential flow solution of tbe 
fluid-elastic force matrix. In it, the diagonal sub-Batrlces [o'] 
and[ ff ] are non-zero and off-diagonal sub-matrices [o" ]and fi"] 
are zero, troro Eqn. (2.14) It is clear that buckling instability 
is controlled by [a"^nd [rjmatrices and dynamic Instability 
by [o |a»l ft" Joiatrloes. Thus as per potential flow theory, 
orily buckXirij’ instability can occur. 


iable 4 gives the potential flow solution of the fluid- 
dampln^’'; force matrix. Diagonal subnaatrices [a']and [ ]are zero 
and off-diagonal sub-matrices [o'']and [t'Jare nonzero in it. 

Prom iiqn, (2,14)^ it is clear that the dynamic instability is 
controlled by [a^Jand [^' ] matrices which are zero as per 
potential jClow theory. 


It clearly indicates that some of the assumptions made in 
calculating the fluid forces by potential flow theory are incorrect. 
As already mentioned, the tube movements are predominantely in 
lift direction for fluid -damping controlled instability. This 
phenomenon is similar to galloping of transmission lines. So the 
vortex shedding behind the tubes may be a source for this 
Instability and the assumption of irrotational flow is incorrect. 

In case of fluid -elastic controlled instability, the 
cross flow velocity is high. At this velocity, separation of 
flow is possible. Hence the ass-umption of no separation of flow 
is Incorrect in finding the fluid forces by potential flow theory. 
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Leaver and V/eaver l[ 27 ] , noted that, when a tube 

moves from its equilibrium position, the stream lines of the 
flow will change but this redistribution will not take 
immediately. There is a phase lag for the distribution of flow 
due to fluid inertia. This poinu should be included in finding the 

fluid forces. 


bo the method of calculating the fluid forces needs a 
lot of improvement, 

0.6 VIBIUTION OF TUBE BANKS SUBJECTED TQ_THE 

'CfiOas plW'AiC) axial interml flw t 

In this section, the effect of axial internal flow, V, 
on the critical reduced cross flow velocity, U^. , is studied for 
a row of 5 tubes only. Ibe fluid forces acting on the tube row 
due to cross flow are obtainable from the potential flow theory 
and experimenbaily measured data. Both of them have been used^ 
for this analysis. Below the critical axial flow velocity, V^, 
only, the effect is studied (critical axial flow velocity, V , 
is one at which first instability occurs for a tube suojeoted t 
axial Interml flow only). Because .above this axial internal 
flow velocity, the internal flow alone can cause me instability. 

First, the row of five uniformly spaced simply -supported 

4 . 4 . 4=1 flow theory is studied. The arrangements 

tubes based on potential flow xn y 

^4 A-n OY 1 of flow are shown in Fig. 9. g^P 
of the tubes and direction of xio 

4 4-4 r/R inside fluid mass ratio, B and damping ratio, 

to radiiis ratio, G/R* insiae xi-ux 

A n 55 and 0.02 respectively. The results o 
taken as 0.66, 0.55 ana 


t, are 
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the critical reduced cross flow velocity, U^;, of the tube 1 
row are given in B'ig. 9 as a function of maL clamping 
parameter, 6^, for VsO, 2.0, 2,5. The critical reduced cross 
flow velocity, U*, decreases with the increase in axial 
internal flow velocity, V. Fig, 10 shows, the critical reduced 
vross flow velocity, as a function of axial internal flow 
velocity, V, for 6^= 1,0, 5,0, 50.0, From these figures, it is 
also noted that the effect on the critical reduced cross flow 
velocity is small for low axial flow velocity, V and it is large 
when it approaches the critical axial flow velocity, V*. The 
fixed— fixed tube— row is also studied and the results are given in 
Fig, 11 and 12, Similar trends are observed in this case also, 

'With the experimentally measured fluid forces, the 
critical reduced cross flow velocities of simply -supported tube 
are also calculated and the resiiits are given in Fig, 13 as a 
function of for V= 0.0, 2.0, 2.5. Fig. 14 shows the critical 
reduced cross flow velocity, U^, as a function of axial internal 
flow velocity for «= 1,0, 5,0, 50.0. From these figures it is 
noted that the axial internal flow velocity reduces the critical 
cross flow velocity in the region of > 15. This can be explained 
as follows-Paidoussis [ 3 ] showed for a simply- supported pipe 

With internal flow, the frequency decreases with axial flow 
velocity and becomes zero at critical velocity (w) and buckles , 
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This shows that stiffness goes on decreasing with increase 
in velocity and damping has no effect. Thus internal flow 
has effect in the fluid -elastic controlled instability region 
>15. Fixed -Fixed tubes are also studied and the results 
are shown in Fig, 15 and 16, Similar trends are also observed 
in this case. 



Critical Reduced Cross Flow Velocity, 



Fig.9 Critical Cross Flow Velocity Ur. as a Function of 6m for 
a Row of 5 Simply- Supported Tubes Subjected to theCc 
mbined Flows( Based on Potential Flow Fluid Forces.! 




Critical Reduced Cross Flow Velocity, Ur. 



Flg.10 Critical Cross Flow Velocity , Ur, os a Function of V 
for a Row of 5 Simply -Supported Tubes Subjected 
to the Combined Flows (Based on Potential Flow 



Critical Reduced Cross Flow Velocity, 


50 



Mass -Damping Parameter, 6m 


Fiatl Critical Cross Flow Velocity Ur. as a Function for a 

Row of 6 Fixed -Fixed Tubes Subjected to the Combined 
Flows ( Based on Potential Flow Fluid Forces.) 




Critical Reduced Cross Flow Velocity, Ur . 


1 1. T -ija 



Flg.12 Critical Cross Flow Velocity U*as a Function of V 
for Q Row of 5 Fixed-Fixed Tubes Subjected to the 
Combined Flows (Based on Potential Row Fluid Fo- 

r 




Critical Reduced Cross Flow Velocity, Ur . 
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FiaC Critical Cross Flow Velocity, Ur. os a Function of 6m tor 
a Row of 5 Simply-Supported Tubes Subjected to the 
Combined Flows (Based on Experimentally Measured 
Fluid Forces) 




Critical Reduced Cross Flow Velocity, 
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Fig. 14 Critical CrossFlow Velocity Ur, as a Function of Vfor 
a Row of 5 Simply -Supported Tubes Subjected to the 
Combined Flows (Based on ExperimentaSy Measured 

m. \ 



* 

Fig.15 Critical Cross Flow Velocity Ur, os a Function of bm tor a 
Row of 5 Fixed -Fixed Tubes Subjected to the CoTnbined 
Flows ( Based on Experimentally Measured Fluid Forces.) 



Critical Reduced Cross Flow Velocity, Ur . 



■X* 

Fig,1 6 Critical Cross Flow Velocity Ur, as a Function of 
V for a Row of 5 Fixed -Fixed Tubes Subjected t> 
the Combined Flows (Based on Experimentatl y 
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CHAPTER-IV 

CONCLUSIONS 


In the present work, the equations of motion of a 
tube bank vibrating in the combined cross and axial internal 
flow have been derived. These equations are transfermed into 
a set of simultaneous ordinary differential equations, using FEM. 
From these equations, results are obtained for tube banks in 

(i) stationary outside fluid, 

(ii) cross flow and 

(iii) combined cross and axial internal flow. 

for some of the parameters like outside fluid mass ratio, gap 
to radius ratio and number of tubes with different configuration. 
Based on these results and discussions presented in Chapter III, 
the following conclusion are arrived at 

1 . The improved mathematical model for analysing 
combined cross and axial internal flow is a general one . This 
model has been used for analysing the tube banks in (a) stationai 
outside fluid, (b) cross flow and (c) combined cross and axial 

internal flow. 

2, The finite element solution for the vibration of 
simply -supported tube-row vibrating in stationary outside fluid 
agrees well that of closed form solution by Chen [ZO], which 
validates the FEM formulation developed in this work. 
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3» In cross flow, the results obtained by the fluid 
forces of potential flow theory are unrealistic . Hence the 
potential flow theory is unacceptable for finding the fluid 
forces in cross flow. Further improvements are needed in finding 
these forces by analytical. 

4, Results obtained for the combined cross and axial 
internal flow of the simply -supported and fixed -fixed tube -row, 
using measured fluid forces, indicate that the axial flov/ 
velocity reduces the critical cross flow' velocity in the fluid- 
elastic controlled instability region (U^ > 15 and 6^ >4,5). 

There is a considerable scope of research being conducted 
in the combined cross and axial internal flow induced vibration 
in particular the mechanism associated with the onset of 
instabilities , 
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